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ABSTRACT 

We measure the two-point correlation function £(r p ,ir) in a sample of 2219 galaxies between z = 
0.7 — 1.35 to a magnitude limit of Rab — 24.1 from the first season of the DEEP2 Galaxy Redshift 
Survey. From £(r p ,7r) we recover the real-space correlation function, £(r), which we find can be 
approximated within the errors by a power-law, £(r) = (r/ro)~ 7 , on scales ~ 0.1 — 10 h^ 1 Mpc. In a 
sample with an effective redshift of z c g = 0.82, for a ACDM cosmology we find ro = 3.53 ± 0.81 hr 1 
Mpc (comoving) and 7 = 1.66 ± 0.12, while in a higher-redshift sample with z e g = 1.14 we find ro 
= 3.12 ± 0.72 hr 1 Mpc and 7 = 1.66 ± 0.12. These errors are estimated from mock galaxy catalogs 
and are dominated by the cosmic variance present in the current data sample. We find that red, 
absorption-dominated, passively-evolving galaxies have a larger clustering scale length, ro, than blue, 
emission-line, actively star- forming galaxies. Intrinsically brighter galaxies also cluster more strongly 
than fainter galaxies at z ~ 1. Our results imply that the DEEP2 galaxies have an effective bias 
b = 0.96 ± 0.13 if as dm = 1 today or b = 1.19 ± 0.16 if as dm = 0.8 today. This bias is lower than 
what is predicted by semi-analytic simulations at z ~ 1, which may be the result of our i?-band target 
selection. We discuss possible evolutionary effects within our survey volume, and we compare our 
results with galaxy clustering studies at other redshifts, noting that our star-forming sample at z ~ 1 
has very similar selection criteria as the Lyman-break galaxies at z ~ 3 and that our red, absorption- 
line sample displays a clustering strength comparable to the expected clustering of the Lyman-break 
galaxy descendants at z ~ 1. Our results demonstrate that galaxy clustering properties as a function 
of color, spectral type and luminosity seen in the local Universe were largely in place by z ~ 1. 

Subject headings: galaxies: statistics, distances and evolution 
universe — surveys 
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1. INTRODUCTION 

Understanding the nature of large-scale structure in 
the Universe is a key component of the field of cosmol- 
ogy and is vital to studies of galaxy formation and evolu- 
tion. The clustering of galaxies reflects the distribution 
of primordial mass fluctuations present in the early Uni- 
verse and their evolution with time and also probes the 
complex physics which governs the creation of galaxies 
in their host dark matter potential wells. 

Since the first redshift surveys, the two-point correla- 
tion function, £(r), has been u sed as a measure of th e 
strength of galaxy clustering l)Davis fc Peebles! 119831. 
£(r) is relatively straightforward to calculate from pair 
counts of galaxies, and it has a simple physical interpre- 
tation as the excess probability of finding a galaxy at a 
separation r from another randomly-chosen galaxy above 
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that for an unclustered distribution (Peebles 1980). Lo- 
cally, £(r) follows a power-law, £( r)= (r/rp)~ 7 , on scales 
~ 1 - 10 h' 1 Mpc with 7 ~ 1-8 (iDavis fc Peebles! Il983l 
Ide La poarent. Geller. fc Huchral 119881 iTuckeretld] 
119971 iZehavi et alJ 120021: lHawkins etafl 12003(1 . The 
scale-length of clustering, r , is the separation at 
which the probability of finding another galaxy is 
twice the random probability. Locally ro is mea- 
sured to be ~ 5.0 hr 1 Mpc for optically-selected 
galaxies but depends strongly on galaxy morphol- 
ogy, color, ty pe and luminos ity llDavis fc Gellerj 1976 
' DresslerllT980l lLovedav et all 119951 iHermit et all 
Willmer. da Costa. & Pellegrini 1998: iNorbere 



199C 

_.. _ . et all 

20011 IZehavi et alJl2002t iMadgwick et al.l l: 

The spatial clustering of galaxies need not trace the 
underlying d istribut i on of dark matter. This was first 
discussed by iKaiserl l)1984l) in an attempt to reconcile 
the different clustering scale lengths of field galaxies and 
rich clusters, which cannot both be unbiased tracers of 
mass. The galaxy bias, b, is a measure of the cluster- 
ing in the galaxy population relative to the clustering 
in the underlying dark matter distribution. It can be 
defined as the square root of the ratio of the two-point 
correlation function of the galaxies relative to the dark 
matter: b = (C/Cdm) 1 ^ 2 , either as a function of r or de- 
fined at a specific scale (see Section 5.1). Observations of 
galaxy clustering have shown that the galaxy bias can be 
a function of morphology, type, color, luminosity, scale 
and redshift. 
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Using galaxy morphologies. lLovedav et all l)1995j) find 
that early-type galaxies in the Stromlo-APM rcdshift 
survey are much more strongly clustered than late-type 
galaxies. Their early- type sample has a larger correla- 
tion length, ro, and a steeper slope than late-type galax- 
ies. However, iWillmer. da Costa, fc Pellegrini (1998) 
show using data from the Southern Sky Redshift Sur- 
vey (SSRS2, da Costa et al. 1998) that in the absence 
of rich clusters early-type galaxies have a relative bias 
of only ~ 1.2 compared with late- type galaxies. In their 
sample, red galaxies with (B — R)q > 1.3 are significantly 
more clustered than blue g alaxies, with a relative bias of 
~ 1.4. iZehavi et al.1 l)2002l) also studied galaxy clustering 
as a function of color, using data from the Sloan Digital 
Sky Survey (SDSS, York et al. 2000 ) Early Data Re- 
lease, and find that red galaxies (u* — r* > 1.8) have a 
larger correlation length, ro, a steeper correlation func- 
tion, and a larger pairwise velocity dispersion than blue 
galaxies. They also find a strong dependence of cluster- 
ing strength on luminosity for magnitudes ranging from 
M* + 1.5 to M* - 1.5. Galaxy clustering for different 
spectral types in the 2dF Galaxy Redsh ift Survey (2dF- 
GRS, C olless et al. 2001 ) is reported bv lMadgwick et all 
( 2003a); absorption- line galaxies are shown to have a rel- 
ative bias ~ 2 times that of emission-line galaxies on 
scales r = 1 hr 1 Mpc, declining to unity on larger scales. 
Absorption-line galaxies have a steeper correlation slope 
and a larger pairwise velocity dispersion. All of these re- 
sults indicate that red, absorption-line, early-type galax- 
ies are found predominantly in the more massive virial- 
ized grou ps and clusters i n whi ch the random velocities 
are large. iNorberg et al.l l)2001j) report that the correla- 
tion length of optically-selected galaxies in the 2dFGRS 
depends weakly on luminosity for galaxies fainter than 
L* , the typical luminosity of a galaxy, but rises steeply 
with luminosity for brighter galaxies, with the most lumi- 
nous galaxies being three times more clustered than L* 
galaxies. These results from local z ~ surveys indicate 
that the strength of galaxy clustering is quite sensitive 
to different galaxy properties. 

A critical test of both cosmological and galaxy evolu- 
tion models is the redshift-dependence of galaxy cluster- 
ing. The evolution of the dark matter two-point correla- 
tion function, £dm (r, t) , can be calculated readily and is 
strongly dependent on cosmology. In high-density mod- 
els the clustering strength grows rapidly, while ACDM 
models show a more gradual evolution (e.g., Jenkins et 
al. 1998, Ma 1999) . However, the evolution of the 
galaxy two-point correlation function, £(r, t), depends on 
the evolution of both the underlying dark matter dis- 
tribution and the galaxy bias, which is expected to in- 
crease with redshift. Applying semi-analytic modelling 
of gala xy formation and evolution to dark matter simu- 
lations, iKauffmann et all <|1999bfl present ACDM mod- 
els with ro ~ 4 hT x Mpc for the galaxy distribution at 
z = 1 compared to ro ~ 5.2 h^ 1 Mpc locally. They pre- 
dict a galaxy bias of b ~ 1.2 at z = 1 for galaxies with 
Mb < —19 + 5 log (h) but also find that the galaxy bias 
can be a strong function of luminosity, s tar-formation 
rate, galaxy type, and sample selection. iBenson et al.l 
( 2001), who also apply semi-analytic modelling to ACDM 
dark matter simulations, predict a bias of b = 1.5 at z — 1 
for galaxies with Mb < —19.5 + 5 log (h). They also 
predict a similar morphology-density relation at z = 1 to 



that seen locally. 

Previous redshift surveys which have attempted to 
probe intermediate redshifts from z = — 1 have been 
hampered by small volumes and the resulting severe cos- 
mic variance. Results from the Canada-France Redshift 
Survey (CFRS, LeFevre et al. 1996) are based on ~ 600 
gala xies covering 0.14 degrees 2 , the Norris Redshift Sur- 
vey l|Small et alll99 9i) spars ely samples 20 degrees 2 with 
a surv ey of ~ 800 galaxies, iHogg. Cohen, & Blandford 
(2000) report on a sample of ~ 1200 galaxies in two 
very small fields, incl uding the Hubble Deep Field, and 
iCarlberg et alJ l)1997[) present a survey of ~ 250 galaxies 
in a total area of 27 arcmin 2 , finding that correlations 
found in their if-band data are generally greater than 
those found by optically-selected surveys. iSmall et alJ 
(1999) compare results from several surveys which have 
measured the correlation length ro in the range z — — 1 
and illustrate well the uncertainties in and discrepan- 
cies between these results. For an open CDM cosmology, 
the estimates of the comoving correlation length vary 
from ~ 2 — 5 h^ 1 Mpc at z ~ 0.4 — 0.6. In particu- 
lar, the CFRS survey found a much smaller correlation 
length at z > 0.4 than the other surveys, which generally 
are consistent with weak evolution between 2=1 and 
z = 0. A sign ificantly larger s urvey was undertaken by 
CNOC2 ijShenherd et all 1200 lfl . who obtained redshifts 
for ~ 5000 galaxies ove r 1.44 degrees 2 . Most relevant 
for our purposes may be lAdelbergerl l)1999j) . who present 
clustering results for a deep R < 25.5, z ~ 1 sample of 
~ 800 galaxies covering a total of 42.5 arcmin 2 in 5 fields; 
they quote a correlation length of ro ~ 3 hr 1 Mpc for 
a ACDM cosmology, implying that their galaxy sample 
is an unbiased tracer of the mass at z ~ 1. However, 
many of these surveys cover very small fields and are 
likely to underestimate the true clustering. There is a 
well-known systematic bias towards underestimation of 
ro in volumes which are small enough that all galaxies 
are part of a single large-scale structure and in which the 
large-scale modes cannot be sampled. Furthermore, cos- 
mic variance will dominate any measure of clustering in 
volumes which a re too small to be representative samples 
of the Universe ijDavis et al.lll985|) . 

Here we present early results on g alaxy clustering i n 
the DEEP2 Galaxy Redshift Survey {Davis et al.l l 2002). 
an R— band selected survey which was designed to study 
the universe at z ~ 1 with a volume and sampling den- 
sity comparable to local surveys. Our intent in this paper 
is to provide an initial measure of the galaxy clustering 
in our survey at using the first season of data and to 
investigate the dependence of the clustering on galaxy 
properties, splitting the sample by color, spectral type, 
and luminosity. To constrain galaxy evolution models, 
we measure the galaxy bias for the sample as a whole and 
the relative bias between subsamples. This is the first of 
several planned papers on galaxy clustering within the 
DEEP2 survey, and here we focus strictly on analysis of 
spatial correlations. Discussion of reds hift-space distor - 
tions will appear in a subsequent paper ((Coil et all2004[) . 
In the data from the first season of observations we mea- 
sured 5042 redshifts with z > 0.6 in three fields with a 
total area of 0.72 degrees 2 . The most complete field cur- 
rently covers 0.32 degrees 2 and includes 2219 galaxies in 
the redshift range z = 0.7 — 1.35, which we use as the 
primary data sample in this paper. 
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The outline of the paper is as follows: in Section 2 we 
briefly describe the survey, provide details of the obser- 
vations, data reduction, and the data sample used here. 
Section 3 outlines the methods used in this paper, while 
Section 4 presents our results, both for the survey sample 
as a whole and for subsamples based on galaxy redshift, 
color, spectral type, and luminosity. In Section 5 we 
discuss galaxy bias and the relative biases between our 
subsamples, and we conclude in Section 6. 

2. DATA 

2.1. The DEEP2 Galaxy Redshift Survey 

The DEEP2 Galaxy Redshift Survey is a three-year 
proje ct using the DEIMOS spectrograph (Fabcr ct al. 

2002) on the 10-m Keck II telescope to survey optical 
galaxies at z ~ 1 in a comoving volume of approximately 
6xl0 6 h~ 3 Mpc 3 . The completed survey will cover 3.5 
degrees 2 of the sky over four widely separated fields to 
limit the impact of cosmic variance. The "1-hour-survey" 
(1HS) portion of the DEEP2 project will use ~ 1 hr ex- 
posure times to measure redshifts for ~ 60, 000 galaxies 
in the redshift range z ~ 0.7—1.5 to a limiting magnitude 
of Rab = 24.1 (all magnitudes in this paper are in the 
AB system; Oke fc Gunn 1983 ). Photometric data were 
taken in B, R and /-bands with the CFH12k camera on 
the 3.6-m Canada-France-Hawaii telescope. Galaxies se- 
lected for spectroscopy must additionally meet a color se- 
lection given approximately by B—R < 2.35(i?— I)— 0.45, 
R — I > 1.15, or B — R < 0.5. This simple color-cut was 
designed to select galaxies at z > 0.7 (details are given 
in Newman et al. 200 4) and has proven e ffective in do- 
ing so. As discussed in lDavis et alJ {2002) this color-cut 
results in a sample with ^90% of the objects at z > 0.7, 
missing only ^5% of the z > 0.7 galaxies. 

Each of the four DEEP2 1HS fields corresponds to a 
volume of comoving dimensions ^ 20 x 80 x 1000 
Mpc in a ACDM model at a redshift of z = 1. To convert 
measured redshifts to comoving distances along the line 
of sight we assume a flat cosmology with Sl m — 0.3 and 
£1a = 0.7. Changing cosmological mod els within the 
range allowed by recent WMAP analysis l)Soergel et al.l 

2003) has only a modest influence on our results. We use 
h — 77o/(100 km s -1 ), and we quote correlation lengths, 
ro, in comoving dimensions of Mpc. 

2.2. Observations and Data Reduction 

This paper uses data from the first observing season 
of the 1HS portion of the DEEP2 survey, from August- 
October 2002. Three of the four DEEP2 fields were ob- 
served with a total of 68 custom-made slitmasks. Each 
mask has on the order of ~ 120 slitlets, with a me- 
dian separation in the spatial direction between targeted 
galaxies of ~ 6", and a minimum of 3". Due to the high 
source density of objects, we are able to obtain spectra 
for ~ 67% of our targets. Three 20-minute exposures 
were taken on the DEIMOS spectrograph with a 1200 
line mm -1 grating for each slitmask, covering a spec- 
tral range of ~ 6400 — 9100 A at an effective resolution 
R ~ 5000. The multiple exposures allow us to robustly 
reject cosmic rays from the data. Many of the slitlets 
in each mask are tilted to align with the major axis of 
the target galaxy to enable internal kinematic studies, 
and as a result we do not dither the telescope between 
exposures. 



The data were reduced using a sophisticated IDL 
pipeline developed at UC-Berkeley, adapted from spec- 
troscopic reduction prog rams developed for the SDSS 
l)Burles fc Schleeell l2004fl . To find the redshift of each 
galaxy, a ^-minimization is used, where the code finds 
minima in \ 2 between the observed spectrum and two 
templates; one is an artificial emission-line spectrum con- 
volved with a broadening function to mimic a 1" slit and 
60 km s _1 internal dispersion. The other template is 
a high signal-to-noise ratio absorption-dominated spec- 
trum which is the average of many thousands of SDSS 
gal axies covering a rest wavelength range o f 2700-9000 
A HEisensteiu et nUfxM iBuTleT fc Schle g elll2T)f)l. The 
five most-likely redshifts are saved and used in a final 
stage where the galaxy redshift is confirmed by human 
interaction. Our overall redshift success rate is >70% 
and displays only minor variation with color and magni- 
tude (< 20%), with the exception of the bluest galaxies 
(R— I < 0.4, B — R ' < 0.5) for which our redshift success 
rate is ~ 35%. These galaxies represent ~ 25% of our 
targeted sample and account for ~ 55% of our redshift 
failures. 

The A3727 A [Oil] doublet redshifts out of our spectral 
range at z ~ 1.44, and it is believed that all of our bluest 
(R — I < 0.4, B — R< 0.5) targeted galaxies for which we 
do not measure a redshift lie beyond this range. These 
galaxies have similar colors and source densities as the 
population at z ~ 2 currently studied by C. Steidel and 
collaborators (private communication) . If these galaxies 
were in our observable redshift window, it is almost cer- 
tain that we would have measured a redshift, given that 
these blue galaxies must have recent star-formation and 
therefore strong emission lines. 

Although the instrumental resolution and photon 
statistics of our data would suggest that we could achieve 
a redshift precision of ~ 10 km s" 1 in the rest frame of 
each galaxy, we find using galaxies observed twice on 
overlapping slitmasks that differences in the position or 
alignment of a galaxy within a slit and internal kinemat- 
ics within a galaxy lead to an effective velocity uncer- 
tainty of ~ 30 km s _1 . 

2.3. Data Sample 

Here we present results from only the most complete 
field, centered at 02 hr 30 min +00 deg, for which we have 
observed 32 slitmasks covering ~ 0.7 degrees by ~ 0.5 
degrees on the sky. We use data only from masks which 
have a redshift success rate of 60% and higher in order 
to avoid systematic effects which may bias our results. 
Figures and |3 show the spatial distribution of galaxies 
on the plane of the sky and the window function for this 
field. The observed slitmasks overlap each other in two 
horizontal rows on the sky. Six of the masks have not 
as yet been observed in this pointing, leading to regions 
with lower completeness. 

While we measure redshifts as high as z — 1.48, for 
this paper we include only galaxies with 0.7 < z < 1.35, 
a range in which our selection function is currently well 
defined. Our sample in this field and range contains 2219 
galaxies, with a median redshift of z = 0.90. At this me- 
dian redshift the typical rest-frame wavelength coverage 
is ~ 3400 - 4800 A. Figure shows the overall redshift 
distribution of galaxies with 0.5 < z < 1.5 in all three 
of our observed fields. There is a rise between redshifts 
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z = 0.7 — 0.8, the result of our probabilistic pre-selection 
of spectroscopic targets expected to have redshifts > 0.7. 
The flux limit of our sample results in the slow decrease 
of the observed objects at higher redshifts; smaller-scale 
variations are due to galaxy clustering. 

In order to compute galaxy correlation statistics, we 
must understand our selection function (p(z), defined as 
the relative probability at each redshift that an object 
will be observed in our sample. In general, the selec- 
tion function can depend on redshift, color, magnitude, 
and other properties of the galaxy population and sur- 
vey selection. Ideally one would compute <j)(z) from the 
luminosity function of galaxies in the survey. For this ini- 
tial study we estimate <f>(z) by smoothing the observed 
redshift histogram of all the galaxies in our sample, tak- 
ing into account the change in volume with redshift. We 
smoothed with a boxcar of width 450 h^ 1 Mpc and then 
used an additional boxcar of width 150 ft -1 Mpc to en- 
sure that there were no residual bumps due to large-scale 
structure. The resulting <fi(z) is shown by the solid line 
in Figure |3| Also shown in this figure are the normalized 
selection functions for the emission-line and absorption- 
line samples discussed later in the paper (see Section 4.4). 
Note that the redshift distribution <j){z) is determined us- 
ing galaxies in all three of our observed fields, not only 
in the field for which we measure £(r p , 7r), which reduces 
effects due to cosmic variance. Use of a preliminary 4>(z) 
constructed from the luminosity function of our sample 
does not change the results presented here. Using mock 
catalogs to test the possible systematic effects due to our 
estimation of 4>(z), we find that the resulting error on ro 
is 5%, signficantly less than that due to cosmic variance. 

3. METHODS 

3.1. Measuring the Two-point Correlation Function 

The two-point correlation function £(r) is a measure of 
the excess probability above Poisson of finding a galaxy 
in a volume element dV at a separation r from another 
randomly chosen galaxy, 



the ILandv fc Szalavl l|1993|) estimator, 



dP = n[l + £(r)]dV, 



(1) 



where n is the mean number density of galaxies. To mea- 
sure £(r) one must first construct a catalog with a ran- 
dom spatial distribution and uniform density of points 
with the same selection criteria as the data, to serve as 
an unclustered distribution with which to compare the 
data. For each data sample we create a random catalog 
with initially > 40 times as many objects with the same 
overall sky coverage as the data and uniform redshift cov- 
erage. This is achieved by applying the window function 
of our data, seen in Figure 1, to the random catalog. 
Our redshift completeness is not entirely uniform across 
the survey; some masks are observed under better con- 
ditions than others and therefore yield a higher success 
rate. This spatially-varying redshift success complete- 
ness is taken into account in the window function. We 
also mask the regions of the random catalog where the 
photometric data had saturated stars and CCD defects. 
Finally, we apply our selection function, <f>(z), so that the 
random catalog has the same overall redshift distribution 
as the data. This results in a final random catalog which 
has > 15 times as many points as the data. 
We measure the two-point correlation function using 
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where DD,DR, and RR are pair counts of galaxies in 
the data-data, data-random, and random-random cata- 
logs, and nu and ur are the mean number densities of 
galaxies in the data and random catalogs. This estimator 
has be en shown to per form as well as the Hamilton esti- 
mator l)Hamiltonl [l993) but is preferred as it is relatively 
insensitive to the size of t he random catalog and han- 
dles edge corrections well ijKerscher. Szanudi. fc Szalav 
l2000l) . 

As we measure the redshift of each galaxy and not 
its distance distortions in £ are introduced parallel to 
the line of sight due to peculiar velocities of galaxies. 
On small scales, random motions in groups and clusters 
cause an elongation in redshift-space maps along the line- 
of-sight known as "fingers of God" . On large scales, co- 
herent infall of galaxies into forming structures causes 
an ap parent contr action of structure along the line-of- 
sight ( Kaiser 1987). While these distortions can be used 
to uncover information about the underlying matter den- 
sity and thermal motions of the galaxies, they complicate 
a measurement of the two-point correlation function in 
real space. Instead, what is measured is £(s), where s is 
the redshift-space separation between a pair of galaxies. 
In order to determine the effects of these redshift-space 
distortions and uncover the real-space clustering proper- 
ties, we measure £ in two dimensions, both perpendicu- 
lar to and along the line of sight. Following iFisher et alJ 
(1994), we define vi and v 2 to be the redshift positions 
of a pair of galaxies, s to be the redshift-space separation 
(vi — V2), and 1 = |(vi + V2) to be the mean distance 
to the pair. We then define the separation between the 
two galaxies across (r p ) and along (w) the line of sight as 



s • 1 

IP 



(3) 



- 7T 2 . (4) 

In applying the ILandv fc Szalavl l)199.'ff ) estimator, we 
therefore compute pair counts over a two-dimensional 
grid of separations to estimate £(r p ,7r). 

In measuring the galaxy clustering, one sums over 
counts of galaxy pairs as a function of separation, nor- 
malizing by the counts of pairs in the random catalog. 
While £(Vp,7r) is not a function of the overall density 
of galaxies in the sample, if the observed density is not 
uniform throughout the sample then a region with higher 
density will contribute more to the total counts of galaxy 
pairs, effectively receiving greater weight in the final cal- 
culation. The magnitude-limit of our survey insures that 
our selection function, 4>(z), is not flat, especially at the 
higher redshift end of our sample, as seen in Figure 2. 
To counteract this, one might weight the galaxy pairs by 
l/(j)(z), though this will add significant noise where 4>{z) 
is low. Wha t is generally used inste ad is the J3 weight- 
ing method l|Davis fc Huchral 1 1982ft . which attempts to 
weight each volume element equally, regardless of red- 
shift, while minimizing the variance on large scales. Us- 
ing this weighting scheme, each galaxy in a pair is given 
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a weight 

Mr) = / Z(s)s 2 ds, (6) 
Jo 

where Z\ is the redshift of the galaxy, t is the redshift- 
space separation between the galaxy and its pair object, 
t = |si — S2I, 4>{z) is the selection function of the sam- 
ple, such that the mean number density of objects in the 
sample is nD<t>{z) for a homogeneous distribution, and J3 
is the volume integral of £(s). We limit r < 20 hr 1 Mpc, 
the maximum r p separation we measure, in order to not 
over-weight the larger scales, which would lead to a nois- 
ier estimate of £(r p , 7r). Note that the weighting depends 
on the integral over £(s), a quantity we want to measure. 
Ideally one would iterate the process of estimating £(s) 
and using the measured parameters in the J3 weighting 
until convergence was reached. Here we use a power- 
law form of £(s) = (s/so)~ 7 , with initial parameters of 
so = 4.4 hr 1 Mpc and 7=1.5. These power-law values 
are in rough accordance with £(s) as measured in our 
full sample. As tests show that the measured £(r p ,7r) i s 
quite insensitive to the assumed values of sq and 7, we do 
not iterate this process. We estimate no to be 0.003/i 3 
Mpc -3 from the observed number density of galaxies in 
our sample in the redshift range z — 0.75 — 0.9. As with 
so and 7, we find that the results are not sensitive to the 
exact value of used. 

3.2. Deriving the Real-Space Correlations 

While £(s) can be directly calculated from pair counts, 
it includes redshift-space distortions and is not as easily 
interpreted as the real-space correlation function, 

which measures only the physical clustering of galaxies, 
independent of any peculiar velocities. To recover £ (r) we 
use a projection of £(r p , 7r) along the r p axis. As redshift- 
space distortions affect only the line-of-sight component 
of £(r p ,7r), integrating over the 7r direction leads to a 
statistic w p (r p ), whic h is independent of red shift-space 
distortions. Following iDavis fc Peebles! l|1983j) . 

/>oo />OC 

w p (r p ) = 2 dir £(r p ,7r) = 2 / dy £(r* +y 2 ) 1, \ 
Jo Jo 

(7) 

where y is the real-space separation along the line of 
sight. If £(r) is modelled as a power-law, £(r) = 
( r / r o) then ro and 7 can be readily extracted from 
the projected correlation function, w p (r p ), using an an- 
alytic solution to Equation [7] 

w P (r p )=r p ^-) r(i) ' (8) 

where T is the usual gamma function. A power-law fit 
to w p (r p ) will then recover ro and 7 for the real-space 
correlation function, £(r). In practice, Equation is 
not integrated to infinite separations. Here we integrate 
to 7r max = 20 hr 1 Mpc, which includes most correlated 
pairs. We use analytic calculations of a broken power- 
law model for £(r) which becomes negative on large scales 
that we are underestimating ro by less than ~2% by not 
integrating to infinity. 



3.3. Systematic Biases due to Slitmask Observations 

When observing with multi-object slitmasks, the spec- 
tra of targets cannot be allowed to overlap on the CCD 
array; therefore, objects that lie near each other in the di- 
rection on the sky that maps to the wavelength direction 
on the CCD cannot be simultaneously observed. This 
will necessarily result in under-sampling the regions with 
the highest density of targets on the plane of the sky. To 
reduce the impact of this bias, adjacent slitmasks are po- 
sitioned approximately a half-mask width apart, giving 
each galaxy two chances to appear on a mask; we also 
use adaptive tiling of the slitmasks to hold constant the 
number of targets per mask. In spite of these steps, the 
probability that a target is selected for spectroscopy is 
diminished by ~ 25% if the distance to its second near- 
est neighbor is less than 10 arcseconds (see Davis et al. 
2004 for details). This introduces a predictable system- 
atic bias which leads to underestimating the correlation 
strength on small scales. 

Some previous surveys have attempted to quantify and 
correct for effects of this sort using the projected corre- 
lation function, w(8), of the sample selected for spec- 
troscopy rel ative to the that of the entire photomet- 
ric sample l|Hawkins et all l2~003|) . Other surveys have 
attempted to correct these effects by giving additional 
weight to observed galaxies which were close to galaxies 
which were not observed or by restricting the sca les on 
which they measure clustering <|Zehavi et al.ll2002fl . It is 
not feasible for us to use measures of w(6), as the line-of- 
sight distance that we sample is large (>1000 hr 1 Mpc) 
and the resulting angular correlations projected through 
this distance are quite small. In addition, the relation 
between the decrease in the 2-dimensional angular cor- 
relations and the 3-dimensional real-space correlations is 
not trivial and depends on both the strength of clustering 
and the redshift distribution of sources. 

In order to measure this bias, we have chosen to use 
mock galaxy catalogs which have similar size, depth, and 
selection function as our survey and which simulate the 
real-space clustering present in our data. We have con- 
structed these mock catalogs f rom the GIF semi-analyti c 
models of galaxy formation of|Kajfffinajm et_aL| (^999a) . 
As described in ICoil. Davis, fc SzapudlpOOll) . we use 
outputs from several epochs to create six mock catalogs 
covering the redshift range z = 0.7 — 1.5. To convert 
the given comoving distance of each object to a redshift 
we assumed a ACDM cosmology; we then constructed a 
flux-limit ed sample which has a similar source density as 
our data. ICoil. Davis. & Szanudil 1)200 If ) presents the se- 
lection function and clustering properties of these mock 
catalogs. 

To quantify the effect of our slitmask target selection 
on our ability to measure the clustering of galaxies, we 
calculate £(r p , w) and w p (r p ) in these mock catalogs, both 
for the full sample of galaxies and for a subsample that 
would have been selected to be observed on slitmasks. 
The projected correlation function, w p (r p ), of objects se- 
lected to be on slitmasks is lower on scales r < 1 h^ 1 
Mpc and higher on scales r > 1 h^ 1 Mpc than the full 
catalog of objects. We find that ro as measured from 
w p (r p ) is overestimated by 1.5% in the targeted sample 
relative to the full sample, while 7 is underestimated by 
4%. Thus our target selection algorithm has a relatively 
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small effect on estimates of the correlation strength that 
is well within the expected uncertainties due to cosmic 
variance. We do not attempt to correct for this effect in 
this paper. 

4. RESULTS 

We show the spatial distribution of galaxies in our 
most complete field with 0.7 < z < 1.35 in Figure 0| 
We have projected through the short axis, correspond- 
ing to declination, and plot the comoving positions of 
the galaxies along and transverse to the line of sight. 
Different symbols show emission-line and absorption-line 
galaxies, classified by their spectral type as discussed in 
Section 4.4. Large-scale clustering can be seen, with 
coherent structures such as walls and filaments of size 
> 20 h^ 1 Mpc running across our sample. There are 
several prominent voids which contain very few galaxies, 
and several overdense regions of strong clustering. The 
visual impression is cons i stent with ACDM c osmologies 
ijKauffmann et alJll999bl IBenson et al.112001^ . An anal- 
ysis of galaxy groups an d clusters in the ea rly DEEP2 
data will be presented bv lOerke et al.1 l)2004j) . 

In this paper we focus on measuring the strength of 
clustering in the galaxy population using the two-point 
correlation function. First we measure the clustering for 
the full sample shown in Figure 0| Given the large depth 
of the sample in redshift, we then address whether it is 
meaningful to find a single measure of the clustering over 
such an extended redshift range, as there may be signifi- 
cant evolution in the clustering strength within our sur- 
vey volume. To investigate evolution within the sample, 
we would like to measure the clustering in limited red- 
shift ranges within the survey; given the current sample 
size, we divide the data into only two redshift subsam- 
ples, studying the front half and back half of the survey 
separately. Finally, we split the full sample by predicted 
restframe (B — R)o color, observed R — I color, spectral 
type, and absolute Mb luminosity, to study galaxy clus- 
tering as a function of these properties at z ~ 1. The 
survey is far from complete, and with the data presented 
here we do not attempt to subdivide the sample further. 
In future papers we will be able to investigate the clus- 
tering properties of galaxies in more detail. 

4.1. Clustering in the Full Sample 

The left side of Figure shows £(r p ,7r) as measured 
for all galaxies in the most complete field of our survey 
in the redshift range z = 0.7 — 1.35. All contour plots 
presented here have been produced from measurements 
of £(r p , 7r) in linear bins of 1 x 1 h^ 1 Mpc, smoothed with 
a 3 x 3 boxcar. We apply this smoothing only for the fig- 
ures; we do not smooth the data before performing any 
calculations. On scales r p < 2 h^ 1 Mpc, the signature 
of small "fingers of God" can be seen as a slight elon- 
gation of the contours in the it direction. Specifically, 
the contours of £ = 2 and £ = 1 (in bold) intersect the 
7r-axis at ~ 3.5 and 5 h^ 1 Mpc, while intersecting the 
r p -axis at ~ 2.5 and 4 hr 1 Mpc, respectively. We leave 
a detailed investig ation of redshift- space distortions to a 
subsequent paper (jCoil et alJl2004|) . 

In order to recover the real-space correlation function, 
£(r), we compute the projected function w p (r p ) by cal- 
culating £(r p ,7r) in log-separation bins in r p and then 
summing over the it direction. The result is shown in the 



top of Figure H3 Errors are calculated from the variance 
of w p {r p ) measured across the six GIF mock catalogs, 
after application of this field's current window function. 
Here w p {r p ) deviates slightly from a perfect power-law, 
showing a small excess on scales r p ~ 1 — 3 h^ 1 Mpc. 
However, the deviations are within the 1-a errors, and 
as there exists significant covariance between the plot- 
ted points, there is no reason to elaborate the fit. From 
w p {r p ) we can compute ro and 7 of £(r) if we assume that 
£(r) is a power-law, using Equation|Hl Fitting w p (r p ) on 
scales r p = 0.1 - 20 hr 1 Mpc, we find r = 3.19 ± 0.51 
and 7 = 1.68 ± 0.07. This fit is shown in Figure[|)]as the 
dotted line and is listed in Tabled 

The errors on ro and 7 are taken from the percentage 
variance of the measured ro and 7 amongst the mock cat- 
alogs, scaled to our observed values. Note that we do not 
use the errors on w p (r p ) as a function of scale shown in 
Figure |H1 which have significant covariance, to estimate 
the errors on tq and 7. The errors quoted on tq and 7 are 
dominated by the cosmic variance present in the current 
data sample. To test that our mock catalogs are inde- 
pendent enough to fully estimate the effects of cosmic 
variance, we also measure the varianc e in new mock cat - 
alogs, which are just being developed lYan et al J |2003). 
These new mocks are made from a simulation with a box 
size of 300 hr 1 Mpc with finely-spaced redshift outputs. 
We find that the error on ro in the new mocks is 12%, 
less than the 16% we find in the GIF mock catalogs used 
here. We therefore believe that the errors quoted here 
fully reflect the effects of cosmic variance. Preliminary 
measurements in two of our other fields are consistent 
with the values found here, within the 1-a errors; in one 
field, with 1372 galaxies, we measure ro = 3.55 hr 1 Mpc 
and 7 = 1.61, while a separate field, with 639 galaxies, 
yields r = 3.22 hr 1 Mpc and 7 = 1.70. 

We have already described above how we use mock 
catalogs to estimate the bias resulting from our slitmask 
target algorithm, which precludes targeting close pairs in 
one direction on the sky. Another method for quantifying 
this effect is to calculate an upper-limit on the clustering 
using a nearest-neighbor redshift correction, where each 
galaxy that was not selected to be observed on a slitmask 
is given the redshift of the nearest galaxy on the plane 
of the sky with a measured redshift. This correction 
will significantly overestimate the correlations on small 
scales, since it assumes that members of all close pairs 
on the sky are at the same redshift, but it should provide 
a strong upper limit on the correlation length ro. Using 
this correction we find an upper limit on ro of 3.78 ±0.60 
h- 1 Mpc and on 7 of 1.80 ± 0.07. 

4.2. Clustering as a Function of Redshift 

In the above analysis, we measured the correlation 
properties of the full sample over the redshift range 
z = 0.7 — 1.35. This is a wide range over which to 
measure a single clustering strength, given both possi- 
ble evolutionary effects in the clustering of galaxies and 
the changing selection function of our survey, as the lumi- 
nosity distance and the rest-frame bandpass of our selec- 
tion criteria change with redshift. In addition to possibly 
'washing out' evolutionary effects within our survey in 
measuring a single clustering strength over this redshift 
range, the changing selection function makes it difficult 
to interpret these results. In this section we attempt to 
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quantify the redshift dependence in our clustering mea- 
surements. 

We begin by estimating the effective redshift of the 
correlation function we have measured. The calculations 
of £( r pi 7r ) presented in the previous section used the J3 
weighting scheme, which attempts to counteract the se- 
lection function of the survey, 4>(z), and give equal weight 
to volumes at all redshifts without adding noise. To cal- 
culate the 'effective' redshift z c g of the pair counts used 
to calculate £(r p ,7r), we compute the mean ^-weighted 
redshift by summing over galaxy pairs: 

eff EiEjjVi Wi(z,T) 2 

where i runs over all galaxies and j over all galaxies 
within a range of separations from i, with r m i„ < t < 
Tmax, where r is the redshift-space separation between 
the pair of galaxies. The weight iWj(z,r) is given by 
Equation El and as both galaxies in the pair are at es- 
sentially the same redshift (to within r max or better) we 
use Wi(z,r) 2 instead of Wi(z,T)Wj(z-, r). Note that this 
effective redshift will depend on the range of separations 
considered. For r m i„ = 1 h Mpc and T max — 2 hr x 
Mpc, z e g = 0.96, while for T min — 14 hr 1 Mpc and 
Tmax = 16 h^ 1 Mpc, z e g = 1.11 for the full sample. 
For this reason, assuming only one effective redshift for 
the correlation function of a deep galaxy sample covering 
a large redshift range cannot accurately reflect the true 
redshift dependence. We do however estimate an ap- 
proximate averaged value for z c g for the galaxy sample 
presented here by summing over all pairs of galaxies with 
r p or 7r <20 h~ x Mpc (r m ,„ = h^ 1 Mpc and r max = 20 
h^ 1 Mpc). This yields z e g = 0.99 for this data sample, 
though we caution that it is not immediately clear how 
meaningful this number is, given the wide redshift range 
of our data. All values of z e g quoted in Table 1 are for 
< r < 20 h~ l Mpc. 

If £(V p ,7r) is calculated without J3 weighting, the raw 
pair counts in the survey are dominated by volumes with 
the highest number density in our sample, namely z = 
0.75-0.9. Without J 3 weighting we find r = 3.67±0.59 
h^ 1 Mpc and 7 = 1.65 ± 0.07. The effective redshift of 
this result is found using Equation [§1 setting w(z) = 1 
for all galaxies, yielding z e g = 0.90. The differences be- 
tween the values of tq and 7 derived with and without J3 
weighting could be the result of evolution within the sur- 
vey sample, cosmic variance, and/or redshift-dependent 
effects of our survey selection. 

It is important to stress that our use of the traditional 
J3 weighting scheme for minimum variance estimates of 
£,(r p ,ir) leads to an effective redshift of the pair counts 
which is a function of the pair separation; the correlations 
of close pairs have a considerably lower effective redshift 
than pairs with large separation. This complicates the 
interpretation of single values of ro and 7 quoted for the 
entire survey. In local studies of galaxy correlations, one 
assumes that evolutionary corrections within the volume 
studied are insignificant and that the best correlation 
estimate will be achieved with equal weighting of each 
volume element, provided shot noise does not dominate. 
J3 weighting is intended to provide equal weight per unit 
volume to the degree permitted by the radial gradient in 
source density, but it complicates interpretation of re- 
sults within a volume for which evolutionary effects are 



expected. It is far better to subdivide the sample volume 
between high and low redshift and separately apply J 3 
weighting within the subvolumes. 

To this end, we divide our sample near its median red- 
shift, creating subsamples containing roughly equal num- 
bers of galaxies with z — 0.7 — 0.9 and z = 0.9 — 1.35. 
The effective redshift for the lower-z sample (averaged 
over all separations) is z e g — 0.82, while for the higher- z 
sample it is z c g = 1.14. The selection function for the 
redshift subsamples is identical to that shown in Figure^ 
cut at z = 0.9. The right side of Figure [5] shows the mea- 
sured £(r p , 7r) for both subsamples. At lower redshifts 
the data exhibit a larger clustering scale length, as might 
be expected from gravitational growth of structure. The 
lower-z sample also displays more prominent effects from 
"fingers of God" . The bottom of Figure shows the re- 
sulting w p (r p ) and power-law fits for each redshift range. 
Note that we fit a power-law on scales r = 0.1 — 20 
Mpc for the higher- z sample but fit on scales r — 0.1 — 6 
h^ 1 Mpc for the lower-z sample, as w p (r p ) decreases sig- 
nificantly on larger scales. We have tested for systematic 
effects which could lead to such a decrease and have not 
found any. With more data we will be able to see if this 
dip persists. The lower-z sample exhibits a larger scale 
length than the higher- z sample, though the difference is 
well within the l-er uncertainties due to cosmic variance 
(see TableHfor details). For each subsample we estimate 
the errors using the variance among the mock catalogs 
over the same redshift range used for the data. 

A positive luminosity-dependence in the galaxy clus- 
tering would lead to an increase in ro measured at larger 
redshifts, where the effective luminosity is greater. How- 
ever, evolutionary effects could offset this effect, if ro 
grows with time. We find no significant difference in 
our measured value of ro for the lower-redshift sample. 
Locally, significant luminosity-dependenc e has been seen 
in th e clustering of data in the 2dFG RS ijNorberg et alJ 
l200l|) and SDSS ijZehavi et al J 12002) and, if present in 
the galaxy population at z — 1, could complicate mea- 
surements of the evolution of clustering within our sur- 
vey volume, given the higher median luminosity of galax- 
ies in our sample at larger redshifts. We investigate the 
luminosity-dependence of clustering in our sample in Sec- 
tion 4.5 and discuss possible evolutionary effects in Sec- 
tion 5.2. 

4.3. Dependence of Clustering on Color 

We now measure the dependence of clustering on 
specific galaxy properties. We begin by creating red 
and blue subsamples based on either restframe (B — 
R)o color or observed R — I color, which is a di- 
rect observable and does not depend on modelling K- 
corre ctions. K-correction s were calculated using a sub- 
set of lKinnev et alJ l)1996j) galaxy spectra convolved with 
the B, R and I CFH12k filters used in the DEEP2 Sur- 
vey. These are used to create a table containing the 
restframe colors and K-corrections as a function of z and 
R — I color. K-corrections are then obtained for each 
galaxy using a parabolic interpolation (for more details 
see Willmer et al. 2004 ). The median K-corrections for 
the sample used here are ~ —0.2 for R(z) — B(0) and 
~ —0.9 for R(z) — R(0), where we apply corrections to 
our observed R— band magnit udes, which are the de epest 
and most robustly calibrated ijNewman et H]|l)04). Af- 
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ter applying these K-corrections, we estimate the galaxy 
restframe {B — R)q colors and divide the sample near 
the median color into red, (B — R)o > 0.7, and blue, 
(B — R)q < 0.7, subsets. We further restrict the sub- 
samples to the redshift range z = 0.7 — 1.25. We fit 
for the selection function, (j)(z), for each subsample sep- 
arately, again using data from all three observed fields 
which match these color selection criteria, and find that 
the resulting selection functions for the red and blue sub- 
samples are similar to each other and similar to that 
shown in Figure |3| For these and all subsamples we use 
J3 weighting in measuring £(r p , 7r)and estimate errors us- 
ing mock catalogs with half the original sample size. In 
this way we attempt to replicate the error due to cos- 
mic variance, which depends on the survey volume, and 
sample size. We found that the error did not increase 
signficantly when using half the galaxies, which indicates 
that cosmic variance is the dominant source of error. 

We find that the red galaxies have a larger correla- 
tion length and stronger "Fingers of God" . This trend is 
not entirely unexpecte d, as previous data a t z ~ 1 have 
show n similar effects l|Carlberg et al.lll997t iFirth et al.l 
2002) , though the volume which we sample is much larger 
and therefore less affected by cosmic variance. The top 
of Figure shows w p (r p ) for each sample; fits to £(r) are 
given in Table 

While restframe colors are more physically meaning- 
ful than observed colors, they are somewhat uncertain 
as K-corrections can become large at our highest red- 
shifts. We therefore divide our full data sample into red 
and blue subsets based on observed R — I color. There 
is a clear bimodality in the distribution of R — I col- 
ors of DEEP2 targets, leading to a natural separation at 
R — I ~ 1.1. However, as there are not enough galaxies 
with R — I > 1.1 to provide a robust result, we instead 
divide the full dataset at R — I = 0.9, which creates sub- 
samples with ~ 4 times as many blue galaxies as red. 
We again construct redshift selection functions, <fi(z), for 
each sample and measure £,(r p ,ir) in the redshift range 
z = 0.7 — 1.25. Again, the redder galaxies show a larger 
correlation length and a steeper slope, though the differ- 
ences are not as pronounced as in the restframe color- 
selected samples; see Tabled 

4.4. Dependence of Clustering on Spectral Type 

We next in vestigate the dependenc e of clustering on 
spectral type. iMadgwick et alJ l)2003b|) have performed a 
principal component analysis (PCA) of each galaxy spec- 
trum in the DEEP2 survey. They distinguish emission- 
line from absorption-line galaxies using the parameter 77, 
the distribution function of which displays a bimodality, 
suggesting a natural s plit in the sample . We us e the same 
division employed bv IMadgwick et all (|2003bj) . who de- 
fine late-type, emission-line galaxies as having 77 > — 13 
and early-type, absorption-line galaxies with 77 < —13. 
Our absorption-line subset includes ~ 400 galaxies in the 
redshift range z = 0.7 — 1.25, while the emission-line sam- 
ple has ~ 4 times as many galaxies. In Figure 0] different 
symbols show the galaxy population divided by spectral 
type. The early-type subset can be seen to reside in the 
more strongly clustered regions of the galaxy distribu- 
tion. The middle of Figure shows w p {r p ) measured 
for our spectral-type subsamples, and best-fit values of 
ro and 7 are listed in Table Absorption-line galax- 



ies have a larger clustering scale length and an increased 
pairwise velocity dispersion. Since 77 correlates well with 
color, this result is not unexpected; the bulk of the early- 
type galaxies have red colors, though there is a long tail 
which extends to the median color of the late- type galax- 
ies. Thus the subsamples based on spectral type are not 
identical to those based on color. The spectral-type is in- 
timately related to the amount of current star formation 
in a galaxy, so that we may conclude that actively star- 
forming galaxies at z ~ 1 are significantly less clustered 
than galaxies that are passively evolving. Interestingly, 
the emission-line galaxies show a steeper slope in the cor- 
relation function than the absorption-line gala xies, which 
is not seen at z ~ l)Madgwick et al.|l200 3a). This will 
be important to investigate further as the survey collects 
more data. 

4.5. Dependence of Clustering on Luminosity 

We also split the full sample by absolute Mb mag- 
nitude, after applying K-corrections, to investigate the 
dependence of galaxy clustering on luminosity. We di- 
vide our dataset near the median absolute magnitude, 
at M B = -19.75+5 log (h). Figure |S] shows the selec- 
tion function for each subsample. Unlike the previous 
subsets, here the selection functions are significantly dif- 
ferent for each set of galaxies; 4>(z) for the brighter ob- 
jects is relatively flat, while <fi(z) for fainter galaxies falls 
steeply with z. The bottom of Figure [7| shows w p (r p ) for 
each subsample. We fit w p (r p ) as a power-law on scales 
r p ~ 0.15 — 4 h' 1 Mpc and find that the more luminous 
galaxies have a larger correlation length. On larger scales 
both samples show a decline in w p (r p ), but the brighter 
sample shows a steeper decline; fits are listed in Tabled 

In this early paper, using a sample of roughly 7% 
the size we expect to have in the completed survey, we 
have restricted ourselves to considering only two subsam- 
ples at a time. As a result, in our luminosity subsam- 
ples we are mixing populations of red, absorption-line 
galaxies, which have very different mass-to-light ratios 
as well as quite different selection functions, with the 
star-forming galaxies that dominate the population at 
higher redshifts. However, the two luminosity subsam- 
ples in our current analysis contain comparable ratios of 
emission-line to absorption-line galaxies, with ~ 75% of 
the galaxies in the each sample having late-type spec- 
tra. In future papers we will be able to investigate the 
luminosity-dependence of clustering in the star-forming 
and absorption-line populations separately. 

5. DISCUSSION 

Having measured the clustering strength using the real- 
space two-point correlation function, £(r), for each of the 
samples described above, we are now in a position to 
measure the galaxy bias, both for the sample as a whole 
at z ~ 1 and for subsamples defined by galaxy prop- 
erties. We first calculate the absolute bias for galaxies 
in our survey and then determine the relative bias be- 
tween various subsamples. Using these results, we can 
constrain models of galaxy evolution and compare our 
results to other studies at higher and lower redshifts. 

5.1. Galaxy Bias 

To measure the galaxy bias in our sample, we use 
the parameter <r 8 NL , defined as the standard deviation 
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of galaxy count fluctuations in a sphere of radius 8 h 
Mpc. We prefer this quantity as a measure of the cluster- 
ing amplitude over using the scale-length of clustering, 
ro, alone, which has significant covariance with 7. We 
can calculate ct 8 nl from a power-law fit to £ (r) using the 
formula, 

^-™ij£s5$' (10) 

where 

72 

J2(7) = (3- 7 )(4- 7 )(6- 7 )27 < U > 

IjPeeblesI 11980ft . Note that here we are not using the 
linear-theory <r 8 that is usually quoted. Instead, we are 
evaluating £(r) on the scale of 8 Mpc in the non- 
linear regime, leading to the notation <r 8 NL . 
We then define the effective galaxy bias as 



°8 DM 

where <r 8 NL is for the galaxies and <t 8 N q M is for the dark 
matter. Our measurements of <J 8 NL for all data samples 
considered are listed in Tabled Errors are derived from 
the standard deviation of <r 8 NL as measured across the 
mock catalogs. 

The evolution of the dark matter clustering can be 
predicted readily using cither N-body simulations or an- 
alytic theory. Here we compute a g N np from the dark 
matter simulations of lYan et al.l l)2003T) at the effective 
redshifts of both our lower-z and higher-z subsamples. 
We use two ACDM simulations where the linear <r 8 dm 
at z = 0, defined by integrating over the linear power 
spectrum, is equal to 1.0 and 0.8. This is the cr 8 which 
is usually quoted in linear theory. For convenience, we 
define the parameter s 8 = er 8 dm( z = 0)- I n both of 
these simulations, we fit £(?")dm as a power-law on scales 
r ~ 1 — 8 hr 1 Mpc, and from this measure <r 8 jj M using 
ea uation 1 1 01 above . For the simulation with s 8 = 1.0, we 
measure ct^dm = 0.70 at z — 0.83 and <t^ m = 0.60 
at z = 1.18, while for the simulation with s 8 = 0.8, we 
measure u 8 d M = 0.56 at z = 0.83 and Cg^, = 0.49 at 

2=1.18. 

Our results imply that for s 8 = 1.0 in a ACDM cos- 
mology the effective bias of galaxies in our sample is 
b = 0.96 ± 0.13, such that the galaxies trace the mass. 
This would suggest that there was little or no evolution 
in the galaxy biasing function from z = 1 to and could 
also imply an early epoch of galaxy formation for these 
galaxies, such that by z ~ 1 they have become rela- 
tively unbiased. However, if s 8 = 0.8 then the effec- 
tive galaxy bias in our sample is b = 1.19 ± 0.16, which 
is more consistent with predict ions from semi-analytic 
models ( Kauffman n et al.lfl9 99b) . Generally, we find the 
net bias of galaxies in our sample to be b ~ l/s 8 . 

The galaxy bias can be a strong function of sample se- 
lection. One explanation for the somewhat low clustering 
amplitude we find may be the color selection of the sur- 
vey. Our flux-limited sample in the i?-band translates 
to bands centered at A = 3600 A and 3100 A at red- 
shifts z =0.8 and 1.1, respectively. The flux of a galaxy 
at these ultraviolet wavelengths is dominated by young 



stars, and therefore our sample could undercount galax- 
ies which have had no recent star formation, while pref- 
erentially selecting galaxies with recent star formation. 
The DEEP 2 sample selection may be similar to IRAS- 
selected low- z galaxy samples in that red, old stellar pop- 
ulations are under-represented (however, our UV-bright 
sample is probably less dusty than the IRAS galaxies). 
IRAS-selected samples are known to have a diminished 
correlation amplitude and undercount dense regions in 
cluster cores (e.g. Moore et al. 1994). We are ac- 
cumulating A"-band imaging within the DEEP2 fields, 
which we can use to study the covariance of AT-selected 
samples with our i?-selected sample in order to gain a 
better understanding of the behavior of £(r) at z ~ 1. 
iCarlberg et alJ l)1997j) find that their AT-selected sample 
at z ~ 0.3 — 1 generally shows stronger clustering than 
optically-selected samples at the same redshifts, and we 
expect the same will hold true for our sample. 

5.2. Evolution of Clustering Within Our Survey 

The DEEP2 survey volume is sufficiently extended in 
the redshift direction that we expect to discern evolu- 
tionary effects from within our sample. For example, the 
look-back time to z = 0.8 is 6.9 Gyrs in a ACDM cosmol- 
ogy (for h = 0.7), but at z = 1.2 the look-back time grows 
to 8.4 Gyrs. As discussed in Section 4.2, measuring the 
clustering strength for the full sample from z = 0.7 — 1.35 
is not entirely meaningful, as there may be significant 
evolutionary effects within the sample, and the results 
are difficult to interpret given the dependence of the ef- 
fective redshift on scale. We therefore divide the sample 
into two redshift ranges and measure the clustering in 
the foreground and background of our survey. However, 
with the data available to date, our results must be con- 
sidered initial; we hope to report on a sample ~ 20 times 
larger in the next few years. Note also that as our sample 
size increases and we are better able to divide our sample 
into narrower redshift ranges, the dependence of z e g on 
scale will become much less important. 

The decreased correlations observed in the higher- 
redshift subset within the DEEP2 sample might be con- 
sidered to be the effect of an inherent diminished clus- 
tering amplitude for galaxies in the more distant half of 
the survey. Indeed, the mass correlations are expected 
to be weaker at earlier times, but we expect galaxy bias- 
ing to be stronger, so that the galaxy clustering may not 
increase with time as the dark matter distribution does. 
As a complication, at higher redshifts we are sampling 
intrinsically brighter galaxies due to the flux limit of the 
survey, and there is a significant dependence of clustering 
strength on luminosity in our data. Our lower-redshift 
sample has an effective luminosity of Mb = —19.7+5 log 
(h), while for the the higher-redshift galaxies the effec- 
tive luminosity is Mb = —20.4+5 log (h). As discussed 
in Section 4.2, this luminosity difference would lead to an 
increase in ro measured for the higher-redshift sample, if 
there was no intrinsic evolution in the galaxy clustering. 
Additionally, at higher redshifts our i?-band selection 
corresponds to even shorter restframe wavelengths, yield- 
ing a sample more strongly biased towards star-forming 
galaxies. 

5.3. Comparison with Higher and Lower Redshift 

Samples 
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Galaxies which form at high-rcdshift are expected to 
be highly-biased tracers of the unde rlying dark mat- 
ter density field ijBardeen et alJ 119861): this bias is ex- 
pected to then decrease with time llNusse r fc Davisil994t 
IMo fc White1IT99?i nWiark fc Peebles! 11 998ft . Tf galax- 
ies are born as rare peaks of bias bo in a Gaussian noise 
field with a preserved number density, their bias will de- 
cline with epoch according to b = + 1, where D is 
the linear growth of density fluctuations in the interval 
since the birth of the objects. This equation shows that if 
galaxies are highly biased tracers when born, they should 
become less biased as the Universe continues to expand 
and further structure forms. This has been the usual 
explanation for the surprisingly large clustering ampli- 
tude reported for Lyman-break galaxies at z ~ 3. They 
have a clustering scale length comparable to optically- 
selected galaxies in the local Universe, but the dark mat- 
ter should be much less clustered at that epoch, imply- 
ing a bias of &t,vR = 4 .0 ± 0.7 for a ACDM cosmology 
ijAdelberger et al.l H9981. The 2dFGRS team has shown 
that the bi as in their &j-selected sample is consist ent with 
&2DF = 1 (|Verde et aijl200l iLahav et all 12002ft . Given 
these observations of b = 4 at z ~ 3 and b = 1 at z ~ 0, 
one might expect an intermediate value of b at z ~ 1, 
assuming that all of these surveys trace similar galaxy 
populations. However, different selection criteria may be 
necessary to trace the same galaxy population over var- 
ious redshifts. 

Our subsample of star-forming, emission-line galaxies 
has similar selection criteria as recent studies of galaxies 
at z ~ 3. The Lyman-break population has been selected 
to have strong UV luminosity and therefore high star- 
formation rates. The spectroscopic limit of the Lyman- 
break sample is R ~ 25.5, which is roughly equivalent 
to R = 23.5 at z ~ 1, while the DEEP2 survey limit is 
R = 24.1, so that roughly similar UV luminosities are 
being probed by these studies. With a sample of ~ 700 
Lyman-break galaxies at z ~ 3 lAdelberger et alJ (|20Q3Tl 
measure a correlation length of ro = 3.96±0.29 hr 1 Mpc 
with a slope of 7 = 1.55±0.15. At z g = 0.99 we measure 
a somewhat lower correlation length of ro = 3.19 ± 0.51 
h^ 1 Mpc and a slightly steeper slope of 7 = 1.68 ± 0.07, 
implying that star-forming galaxies at z ~ 1 are not as 
strongly biased at those at higher redshifts. The slope of 
the correlation function is expected to increase with time, 
as seen here, as the underlying dark matter continues to 
cluster, resulting in more of the mass being concentrated 
on smaller scales. In constraining galaxy evolution mod- 
els, however, it is important to note that while these are 
measures of similar, star-forming populations of galax- 
ies at z ~ 3 and z ~ 1, the Lyman-break galaxies are 
not progenitors of the star-forming galaxies at z ~ 1. 
Using the linear approximations of iTegmark fc Peebles! 
(1998) one would expect the Lyman-break galaxies to 
have a correlation length of ro ~ 5 hr 1 Mpc at z ~ 1 
ljAdelbergerlll999jl . so that the objects carrying the bulk 
of the star formation at z ~ 1 and z ~ 3 are not the 
same. Our population of red, absorption-line galaxies 
have a correlation length of ro ~ 5 — 6 Mpc, similar 
to that expected for the descendants of the Lyman-break 
population at z ~ 1. 

Using recent studies from both 2dF and SDSS, we can 
also compare our results to z ~ surveys. The two- 



point correlation function is relatively well-fit by a power- 
law in all three of these surveys on scales r = 1 — 10 
hr 1 Mpc. SDSS find a correlation length of rp = 6.1 ± 
0.02 h~ x Mpc in their r*-selected sample ijZehavi et alJ 
2002), while 2dF find r = 5.05 ± 0.26 h~ l Mpc in their 
&j-selected survey. These values are significantly larger 
than our measured ro at z ~ 1, in our i?-selected survey. 
The slope of the two-point correlation function may be 
marginally steeper at low redshifts, with 2dF finding a 
value of 7 = 1.67 ± 0.03 and SDSS fitting for 7 = 1.75 ± 
0.03, compared with our values of 7 = 1.66 ± 0.12 at 
Zeff = 0.82 and z e s = 1.14. 

While the highly-biased, star-forming galaxies seen at 
z ~ 3 appear to have formed in the most massive dark 
matt er halos present at that epoch l)Mo. Mao, fc White 
1999) and evolved into the red, clustered population seen 
at z ~ 1, the star-forming galaxies seen at z ~ 1 are not 
likely to be significantly more clustered in the present 
Universe. These galaxies are not highly-biased, and as 
their clustering properties do not imply that they reside 
in proto-cluster cores, they cannot become cluster mem- 
bers at z = in significant numbers. 

5.4. Relative Bias of Subsamples 

Having measured the absolute galaxy bias in our sam- 
ple as a whole, which is largely determined by the details 
of our sample selection, we now turn to relative trends 
seen within our data, which should be more universal. 
Using the various subsamples of our data defined above 
we can quantify the dependence of galaxy bias on color, 
type, and luminosity, and we compare our findings with 
other results at z = — 1. 

We define the relative bias between two samples as the 
ratio of their a s NL : 



As the subsamples are taken from the same volume and 
have similar selection functions, there is negligable cos- 
mic variance in the ratio of the clustering strengths, and 
therefore the error in the relative bias is lower than the 
error on the values of cr g NL individually. To estimate the 
error on the relative bias, we use the variance among the 
mock catalogs (neglecting cosmic variance) which leads 
to a 4% error, and include an additional error of 6% 
due to uncertainties in the selection function, added in 
quadrature. 

We find in the restframe (B — R)o red and blue 
subsamples that b( B -R) >o.7/b(B-R) <o.7 = 1.41 ± 
0.10. This value is quite similar to the relative bi- 
ases seen in local z = samples. In the SSRS2 data 
iWillmer. da Costa, fc Pellegrinil l)1998f) find red galaxies 
with (B — R)o > 1.3 have a relative bias of ~ 1 .4 com- 
pared to blue galaxies, while IZehavi et all l)2002fl report 
that in the SDSS Early Data Release red galaxies (based 
on a split at u* — r* = 1.8) have a relative bias of ~ 1.6 
compared to blue galaxies. We find a a similar value of 
the relative bias at z ~ 1 in our red and blue subsamples, 
implying that a color-density relation is in place at these 
higher redshifts. The observed-frame R — I subsamples 
have a relative bias of &i?-/>o.9/^i?-/<o.9 = 1-29 ± 0.09. 
This value is slightly lower than that of the restframe 
color-selected subsamples, as expected since the observed 
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R — I color of galaxies has a strong redshift-dependence 
over the redshift range we cover, z = 0.7 — 1.25, and is 
therefore less effective at distinguishing intrinsically dif- 
ferent samples. 

Using the PCA spectral analysis we find that the 
absorption-line sample has a clustering length, tq, ~ 2 
times larger than the emission-line sample, with a rel- 
ative bias of fcahsnrntinri — line/remission— lino = 1.77 + 0.12. 

Madgwi ck et alJ lj2003a|) find using 2dFGRS data that lo- 
cally, absorption-line galaxies have a relative bias about 
twice that of emission-line galaxies on scales of r ~ 1 
h^ 1 Mpc, but that the relative bias decreases to unity on 
scales >10 h^ 1 Mpc. The relative bias integrated over 
scales up to 8 /i -1 Mpc is 1.45 ± 0.14 at z ~ 0, similar to 
our result at z ~ 1. Our current data sample is not suf- 
ficiently large to robustly measure the scale-dependence 
of the galaxy bias, though th is should readily be measur- 
able f rom the final dataset. iHogg. Cohen, fc Blandford 
( 2000) find in their survey (with z me a ~ 0.5) that galax- 
ies with absorption-line spectra show much stronger clus- 
tering at small separations, though their absorption-line 
sampl e size is small, with 121 galaxies. iCarlberg et al.l 
( 1997) also report that in the redshift interval z = 
0.3 — 0.9, galaxies with red colors have a correlation 
length 2.7 times greater than bluer galaxies with strong 
[Oil] emission. 

Recently, there have been several studies which have 
found very large clustering strengths for extremely-red 
objects (EROs, R — K > 5) at z ~ 1. Us ing the angular 
correlation function, iDaddi et alJ l)2001l) find a correla- 
tion l ength of r n = 12 ± 3 h' 1 Mpc for EROs at z ~ 1.2, 
while iFirth et al.l l|2002f) find that the correlation length 
isr ~ 7.5 — 10.5/i _1 Mpc. These samples are of rare ob- 
jects which have extreme colors and are quite luminous; 
Firt h et alJ £2002) estimate that their sample is ~ 1 — 1.5 
magnitudes brighter than M*. We find a correlation 
length of ro = 6.61 ± 1.12 for our absorption-line sample, 
which has an effective magnitude of Mb — —20.5 + 5 
\og(h). Given the relatively large clustering strength of 
the absorption-line galaxies in our sample and the lumi- 
nosity difference between our sample and the ERO stud- 
ies, it is possible that in our absorption-line sample we 
are seeing a somewhat less extreme population which is 
related to the EROs seen at z ~ 1. 

We find that the relative bias between luminosity sub- 
samples is &m b <-i9.75/&Mb>-19.75 = 1.24 ±0.14. These 
datasets have significantly different selection functions, 
unlike the previous samples, and the error on the rela- 
tive bias due to differences in cosmic variance between 
the two samples results in an additional 8% error, added 
in quadrature. This is calculated using numerical exper- 
iments utilizi ng the cosmic variance in redshift bins cal- 
culated in lNewman fc Davisl l)2002|) . The bright sample 
has a median absolute magnitude of Mb = —20.3+5 log 
(h), while the faint sample has a median Mb = —19.1+5 
log (/i). As noted in Section 4.5, our luminosity subsam- 
ples include both star-forming galaxies as well as older, 
absorption-line galaxies and cover a wide range in red- 
shift (z — 0.7 — 1.25), possibly complicating interpreta- 
tion of these results. However, both samples have the 
same ratio of early-type to late-type spectra. 



The DEEP2 Galaxy Redshift Survey is designed to 
study the evolution of the Universe from the epoch 
z ~ 1.5 to the present by compiling an unprecedented 
dataset with the DEIMOS spectrograph. With the fi- 
nal sample we hope to achieve a statistical precision of 
large-scale structure studies at z ~ 1 that is compara- 
ble to previous generations of local surveys such as the 
Las Campanas Redshift Survey (LCRS, Shectman et al. 
1996). As we complete the survey, our team will explore 
the evolution of the properties of galaxies as well as the 
evolution of their clustering statistics. 

The correlation analysis reported here is far more ro- 
bust than earlier studies at z ~ 1 because of our greatly 
increased sample size and survey volume. We find val- 
ues of the clustering scale-length, r = 3.53 ± 0.81 h^ 1 
Mpc at z off = 0.82 and r = 3.12 ± 0.72 h~ Y Mpc 
at z c ff = 1.14, which are within t he wide range of 
clustering amplitudes reported ea rlier ijSmall et al.lll999t 
IHogg. Cohen, fc Blandfordl 120001) . This implies a value 
of the galaxy bias for our sample, b — 0.96 ± 0.13 if 
cs dm = 1 today or b = 1.19 + 0.16 if erg dm = 0.8 today, 
which is lower than what is predicted by semi-analytic 
simulations of z ~ 1. Our errors are estimated using 
mock catalogs and are dominated by sample variance, 
given the current volume of our dataset. 

We find no evidence for significant evolution of ro 
within our sample, though intrinsic evolutionary effects 
could be masked by luminosity differences in our red- 
shift subsamples. We see a significantly-increased cor- 
relation strength for subsets of galaxies with red colors, 
early-type spectra, and higher luminosity relative to the 
overall population, similar to the behavior observed in 
low-redshift catalogs. Galaxies with little on-going star 
formation cluster much more strongly than actively star- 
forming galaxies in our sample. These clustering results 
as a function of color, spectral type and luminosity are 
consistent with the trends seen in the semi-analytic sim- 
ulations of iKauffmann et al.l ((1999b) at z = 1, and in- 
dicate that galaxy clustering properties as a function of 
color, type, and luminosity at z ~ 1 are generally not 
very different from what is seen at z = 0. 

The overall amplitude of the galaxy clustering observed 
within the DEEP2 survey implies that this is not a 
strongly biased sample of galaxies. For s 8 = 1-0 (de- 
fined as the linear as dm at z — 0), the galaxy bias is 
b = 0.96 + 0.13, while for s 8 = 0.8, the bias of the DEEP2 
galaxies is b — 1.19 + 0.16. This low bias may result from 
the i?-band selection of the survey, which roughly corre- 
sponds to a restframe /7-band selected sample; the more 
clustered, old galaxies with red stellar populations are 
likely to be under-represented as our sample preferen- 
tially contains galaxies with recent star-formation activ- 
ity. However, the same selection bias applies to Lyman- 
break galaxies studied at z ~ 3, which are seen to be 
significantly more biased than our sample at z ~ 1. 

We are undertaking studies with A'-band data in 
our fields, which should lead to clarification of these 
questions. More precise determinations of the evolu- 
tion of clustering within our survey and the luminosity- 
dependence of the galaxy bias at z ~ 1 awaits enlarged 
data samples, on which we will report in due course. 



6. CONCLUSIONS 
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Table 1. Power-law fits of £(r) for various data samples. a 



Sample 


no. of 


2 range 


2„ff b 


Tn 
' U 


'V 


r range 


NL 
°8 




galaxies 








(h ~ ^ Mpc) 




(h 


Mpc) 




full sample 


2219 


0.7- 


1.35 


0.99 


3.19 ± 0.51 


1.68 ± 0.07 


0.1 - 


- 20 


0.60 ±0.08 


lower 2 sample 


1087 


0.7- 


0.9 


0.82 


3.53 ±0.81 


1.66 ± 0.12 


0.1 


- 6 


0.66 ±0.12 


higher z sample 


1132 


0.9 - 


1.35 


1.14 


3.12 ±0.72 


1.66 ±0.12 


0.1 - 


- 20 


0.59 ±0.11 


(B - R) > 0.7 


855 


0.7- 


1.25 


0.96 


4.32 ± 0.73 


1.84 ±0.07 


0.25 


- 10 


0.79 ±0.12 


(B - R% < 0.7 


964 


0.7- 


1.25 


0.93 


2.81 ±0.48 


1.52 ±0.06 


0.25 - 


- 10 


0.56 ±0.08 


R- I > 0.9 


442 


0.7- 


1.25 


0.90 


3.97 ± 0.67 


1.68 ±0.07 


0.25 


- 8 


0.72 ±0.11 


R-K 0.9 


1561 


0.7- 


1.25 


0.95 


2.89 ± 0.49 


1.63 ± 0.07 


0.25 


- 8 


0.56 ±0.08 


absorption- line 


395 


0.7- 


1.25 


0.86 


6.61 ± 1.12 


1.48 ±0.06 


0.25 


- 8 


1.06 ±0.16 


emission-line 


1605 


0.7 - 


1.25 


0.97 


3.17 ±0.54 


1.68 ±0.07 


0.25 


- 8 


0.60 ±0.09 


M B < -19.75 


899 


0.7- 


1.25 


0.99 


3.70 ± 0.63 


1.60 ±0.06 


0.15 


-4 


0.68 ±0.10 


M B > -19.75 


1088 


0.7- 


1.25 


0.89 


2.80 ±0.48 


1.54 ±0.06 


0.15 


- 8 


0.55 ±0.08 



a These fits have not been corrected for the small bias we find in our mock catalogs due to our slitmask target selection algorithm (see 
Section 3.3 for details). 
b see Equation 9 
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Fig. 1. — Spatial distribution of the full DEEP2 sample of 2219 galaxies projected on to the plane of the sky. 




Fig. 2. — Window function of spectroscopic coverage in our most complete pointing to date. We include the 32 slitmasks which have a 
rcdshift completeness > 60% in our analysis. The greyscale ranges from (white) to 0.86 (black) and corresponds to the probability that 
a galaxy meeting our selection criteria at that position in the sky was targeted for spectroscopy. The total length of this field is 2 degrees; 
only the first ~ 0.7 degrees have been covered thus far. 
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Fig. 3. — Redshift distribution of ~ 5000 galaxies observed in the first season of the DEEP2 survey, covering three separate fields for 
a total of 0.72 degrees 2 . The solid line is a smoothed fit which we use to estimate our selection function, <f>(z), in the redshift range 
0.7 < z < 1.35. The dotted and dashed lines show the normalized selection functions for the emission-line and absorption-line samples, 
respectively. 
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Fig. 4. — Redshift-space distribution of galaxies in early DEEP2 data in our most complete field shown as a function of redshift and 
comoving distance along and projected distance across the line of sight, assuming a ACDM cosmology. We have split the sample by PCA 
classification: black, plus-signs show emission-line galaxies while red, diamond symbols show absorption-dominated galaxies. It is apparent 
that galaxies with early-type spectra are more strongly clustered. 
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Fig. 5. — Left: Contours of the two-dimensional correlation function, £(r p , w), smoothed with a 3 X 3 boxcar, measured for 2219 galaxies 
in the redshift range 0.7 < z < 1.35 in our most complete field to date. The smoothing has been applied only for the figures; it is not used 
in calculations. Contours levels are 0.0 (dashed), 0.25, 0.5, 0.75, 1.0 (bold), 2.0 and 5.0. Right: Contours of £(r p ,ir), smoothed with a 3 X 3 
boxcar, measured for lower-redshift galaxies in our sample (solid contours) and for higher-redshift objects (dashed contours). Contours 
levels are 0.25, 0.5, 1.0 (bold), 2.0 and 5.0. 





Fig. 6. — The projected correlation function, w p (r p ), for the full redshift range (top) and two redshift sub-samples (bottom). The dotted 
lines show power-law fits used to recover ro and 7 of £(r) for each sample, as listed in Table 1. Error bars are computed from the variance 
across mock catalogs and are estimates of the cosmic variance. 



16 



Coil et al. 





Fig. 7. — Top: w p (r p ) measured for red (dashed line) and blue (solid line) subsamplcs, divided according to restframe (B — R)q color. 
The dotted lines show power-law fits used to estimate ro and 7 (see Table f). Middle: w p (r p ) measured for emission-line (rj > —13, solid 
line) and absorption-dominated (rj < —13, dashed line) subsamples classified using PCA. Bottom: w p (r p ) measured for subsamples divided 
according to absolute magnitude, Mb, assuming h = 1. 



DEEP2 Galaxy Clustering 



17 




Fig. 8. — Rcdshift histograms and the heavily smoothed curves used to estimate the selection functions, (fi(z), for subsamples divided 
according to absolute magnitude, Mb, assuming h = 1. 



